Dark Energy and Matter Evolution from Lensing Tomography 
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Reconstructed from lensing tomography, the evolution of the dark matter density field in the 
well-understood linear regime can provide model-independent constraints on the growth function of 
structure and the evolution of the dark energy density. We examine this potential in the context 
that high-redshift cosmology has in the future been fixed by CMB measurements. We construct 
sharp tests for the existence of multiple dark matter components or a dark energy component that 
is not a cosmological constant. These functional constraints can be transformed into physically 
motivated model parameters. From the growth function, the fraction of the dark matter in a 
smooth component, such as a light neutrino, may be constrained to a statistical precision of cr(/) « 
by a survey covering a fraction of sky /sky with redshift resolution Az = 0.1. For 
the dark energy, a parameterization in terms of the present energy density Que, equation of state 
w and its redshift derivative w' , the constraints correspond to a{w) = 0.016/^."^^^^ and a mildly 
degenerate combination of the other two parameters. For a fixed floE, cr(w') = 0.046/^^^^^; for Que 
marginalized cj{w') = Q.GWf'^^J^ . 



I. INTRODUCTION 

The weak gravitational lensing of faint galaxies pro- 
vides the most direct probe of mass distribution in the 
universe (e.g. Q). Moreover the evolution of clustering in 
the mass distribution is arguably the best theoretically- 
grounded probe of the dark energy and dark matter [pl. 
Although observations of weak lensing on large scales M| 
are still in the discovery phase , future wide-field sur- 
veys have the potential to rival the statistical precision 
and cosmological utility of luminosity distance measures 
from supernova surveys Even in the context of 

a precisely-determined homogeneous cosmology, lensing 
measurements are unique in that they probe the cluster- 
ing properties of the dark matter and energy. These are 
fixed by the homogeneous cosmology only under partic- 
ular assumptions of the particle constituents (e.g. cold 
dark matter and scalar field dark energy) j8[ || . 

Much of the critical cosmological information lies in 
the temporal or radial direction. A potential obstacle for 
weak lensing is that the observables are inherently two- 
dimensional. All of the matter along the line-of-sight to 
a distant source contributes to the lensing. For a family 
of cosmological models that is described by a handful of 
parameters, this is not a serious drawback. Lack of radial 
information is largely compensated by a large angular 
dynamic range and external cosmological information. 

Given the lack of compelling models for the dark en- 
ergy and controversies surrounding the phenomenology of 
the dark matter on small scales, it is interesting to con- 
sider a more model-independent approach. Indeed recent 
studies of alternate parameterizations of the dark energy 
have revealed potential ambiguities in the interpretation 
of luminosity distance measurements |ll[ p^. 
To address these issues with weak lensing, recovery of the 
temporal dimension becomes critical. 

With future surveys that possess source photometric 
redshift information, recovery of the lost information is 



possible in principle through tomography. Photometric 
redshift techniques are already being applied and tested 
on current lensing data The full two-point statisti- 
cal information can be regained by cross-correlating the 
lensing observables on all source redshift planes jl^ . This 
method utilizes both the angular clustering and the tem- 
poral evolution of the density field but obscures the na- 
ture and hence the model-dependence of the information. 
Additionally, the joint observables are survey dependent 
and computationally cumbersome to analyze. 

In this paper, we instead isolate the temporal infor- 
mation by applying recently developed techniques to re- 
construct the radial density field itself We 
will further focus solely on the linear regime where pre- 
dictions are well-understood. Even utilizing only this 
theoretically-clean subset of information in the data, fu- 
ture surveys can potentially provide interesting model- 
independent constraints on the properties of the dark 
energy and matter. 

The outline of the paper is as follows. In f^, we discuss 
the me thod for reconstruction and statistical forecasts. 
In j plll , we study constraints on the growth function for a 
fixed homogeneous cosmology and in f: IV the dark energy 
density evolution assuming pure cold dark matter. We 
discuss these results in §0. 



II. TOMOGRAPHIC RECONSTRUCTION 



We begin in §11 A by briefly reviewing the tomographic 
reconstruction of the dark matter density field in a fixed 
background cosmology as studied in We then gen- 
eralize to the case where the cosmology and th e den sity 
field must be jointly recovered from the data in §[IB| and 



revi ew F isher techniques for statistical forecasts in § 11 C . 

we outline the fiducial cosmology and survey 
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parameters used for illustrative purposes in the following 
sections. 
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A. Known Homogeneous Cosmology 



B. Unknown Homogeneous Cosmology 



We consider the data to be the lensing convergence k 
in an angular pixel discretized into bins of source redshift 
composed into a data vector d^. Weak lensing dictates 
that the data is a linear projection of an underlying den- 
sity field plus noise in the convergence measurement n,^ 



with 
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where D is the comoving distance in a flat universe 



Diz) 



dz' 
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with H'^ = 87rG/9tot/3 defining the Hubble parameter and 
5Dj is the width of bin j. Here and throughout sub- 
scripts on matrices are labels, whereas matrix elements 
are denoted as Here A = {Sp/p)/a is the density 

fluctuation in the bin with the growth rate in a matter- 
dominated universe scaled out. Note that these distances 
depend on the assumed cosmology. We assume for now 
that D{z) has been fixed by other observations, e.g. fu- 
ture supernovae surveys, but relax this assumption in the 
following section. 

The minimum variance estimator of the underlying 
density field is given by [fsl 



where the reconstruction matrix 



Here 



Ra« = NaaPIaN-i 



Naa = [P*,aN««P«a]" 



(4) 
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is the noise covariance of the estimator. Note that 
R-AkPkA = I so that the estimator is unbiased. 

The statistical properties of the recovered density field 
contain cosmological information. The recovered density 
field is an average of the density fluctuation over a win- 
dow (or mask) Wi(x) defined by the angular pixel and 
redshift binning. The signal covariance of these density 
averages 
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M/,(k)W^;(k)F(fc), 



where P{k) is the linear power spectrum today, 0i = 
DgTow{zi)/ai is the linear decay rate of the potential field, 
with -Dgrow the linear growth rate of the density field, nor- 
malized so that 02 = 1 in the matter dominated regime, 
and W^i(k) are the Fourier transforms of the windows. 
The two-point statistical properties of the reconstruction 
therefore contain information on the growth rate and un- 
derlying power spectrum of the density field. 



The inversion of Eqn. (^) requires an assumption of 
a distance-redshift relation in PkA- If this relation is 
not fixed by external constraints, then the reconstructed 
density field will be a biased measure of the true density 
field. In this case, both the distance-redshift and growth 
rate must be fit to the data. 

If the assumed D{z) is close to the true D{z), then 
the reconstruction of the previous section still serves as a 
useful representation of the data. The reconstruction ma- 
trix Rak employs a slightly incorrect assumption of the 
projection matrix so that RakPkA I- The statistical 
properties of the estimated density field are encapsulated 
in the noise matrix (^, which remains unchanged, and a 
new signal covariance matrix 
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RakPkA SaA PkA^^Ak • 



(7) 



The two-point statistics now also contain information 
about the distance-redshift relation Diz). 



C. Statistical Forecasts 

The two-point statistical information in the recon- 
struction can be exposed through the familiar Fisher ap- 
proach (e.g. [^). If the parameters that underly the 
two-point correlation are given by a vector Sp, the Fisher 
matrix is given by 



pix 



with the covariance matrix 



tr[C^^CAA,iC^^CAA,j] 
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(8) 



(9) 



Here we have assumed that the convergence is measured 
in A^pix independent pixels. We will often write this factor 
as the total sky coverage in independent pixels iVpix = 
47r/sky/^pix where Api^ is the pixel area in steradians. 
A complete treatment would track the small correlations 
between neighboring pixels on contiguous patches of sky 

@- . . . , . 

The inverse of the Fisher matrix F^^ gives an estimate 

of the covariance matrix Cpp of the measured parameters 
Sp. Note that under a re-parameterization of the space, 
the Fisher matrix transforms as a covariant tensor 



■ pp 



'^pp^pp'^ppi 
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We will use this fact to go from model-independent pa- 
rameterizations of the underlying functions to model- 
dependent ones. 



D. Fiducial Model and Survey 

We take as a fiducial cosmology a flat Otot = 1 universe 
with = 0.3 in cold dark matter, $1;, = 0.05 in baryons. 



f^DE = 0.65 in dark energy; an equation of state of the 
dark energy of w{z) = —1 corresponding to a cosmologi- 
cal constant; a dimensionless Hubble constant, h = 0.65, 
scalar spectral index n — I, and amplitude of the initial 
curvature power spectrum — 4.8 x 10~^ (cts = 0.92, 
see for specific definitions of parameters). Since to- 
mographic reconstruction will mainly be useful for next- 
generation surveys, it is reasonable to assume that CMB 
experiments will by that time have determined many of 
the underlying cosmological parameters to high accuracy. 
For simplicity we will here assume that flch^ , ^th^ j n-j 
and r^tot are completely fixed to their fiducial values. We 
will return to this point in j |v] . These parameters fix the 
shape and high redshift normalization of the potential 
power spectrum and so we will work in the context that 
only the growth function and distance-redshift relation 
need be determined through lensing. 

For definiteness, we will take the fiducial survey to be 
defined with circular pixels of area 1 deg^ so that the 
reconstructed density field is in the linear regime [^ . 
Larger pixels would reduce the number of independent 
pixels and hence increase the sample variance in the sig- 
nal dominated regime. In the figures that follow, we take 
/sky = 0.1 (~ 4000 deg^) but note that all lensing er- 
rors may be rescaled as (/sky/0. For the redshift 
binning, we take Az = 0.1 out to z = 3. With this bin- 
ning, the signal covariance matrix Saa is nearly diagonal 
and the statistics reduce to the evolution of the density 
variance in bins. 

We will assume a convergence noise spectrum of the 
form 

N«« = diag[-ff„jN,] , (11) 

as appropriate for random intrinsic galaxy ellipticities. 
We take n = 3.6 x 10^ gal. deg~^ and 7i.,„s = 0.3 as an 
estimate of the usable galaxies and the shear noise per 
galaxy measured from a space-based platform (A. Re- 
fregier, private communication) and form the number of 
galaxies per bin Ni from a redshift distribution pT| 

^ Dexp[^{D/D,)^], (12) 
az dz 

where D, is set to reproduce a median redshift ^mcd = 1- 
These fiducial survey specifications are chosen to repre- 
sent the upper range of the capabilities of surveys in the 
foreseeable future. 



III. GROWTH RATE 

For a fixed distance-redshift relation and high redshift 
power spectrum, the remaining degrees of freedom in the 
two-point statistical properties of the reconstructed den- 
sity field are contained in the growth function ^(z). To 
study how lensing constrains this function, we begin with 
a model-independent approach. 
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FIG. 1: Principle components of the growth function: (upper) first 
five eigenfunctions of the growth function (thick to thin) plotted 
in discrete redshift bins; (lower) rank ordered inverse rms noise 
A^^^^ (thin) and signal-to-noise ratio. Here and in the following 
figures the fiducial model and survey are assumed throughout with 
/sky = 0-1 S'lid Az = 0.1 redshift bins. 

Consider the set of binned growth rates as the pa- 
rameters to be estimated. With fine-binning of the den- 
sity field Az = 0.1, this leads to estimates that have large 
correlated errors since the reconstruction effectively takes 
differences of noisy data. The long-time scale evolution of 
the density field is faithfully preserved in the reconstruc- 
tion |]l8| . Since structure grows in linear theory on the 
expansion time scale in gravitational instability models, 
this information is sufficient to constrain cosmology. 

To better understand the information contained 
therein, consider the principle component or eigenvec- 
tor decomposition of the Fisher matrix F^^ = SAS* and 
the linear combinations of the data they define 

SA = S*s^, Caa = A-i. (13) 

In other words, the eigenvectors are the redshift repre- 
sentation of a new basis that is complete and yields un- 
correlated, orthogonal measurements with variance given 
by the inverse eigenvalue l/Xi- The largest eigenvalues 
correspond to the minimum variance directions and are 
shown in the upper panel of Fig. |[ The first mode has a 
broad single peak between < z < Zmed = 1 correspond- 
ing to the bell-shaped weight in the lensing projection of 
Eqn. (H) . The higher modes exhibit oscillatory structure 
and capture information on the low order derivatives of 
the growth function around this intermediate redshift as 
well as the region z > z^^cd- The spectrum of eigen- 
values A~^^^, scaled to represent inverse rms noise and 
normalized for /sky = 0.1 are shown in the bottom panel. 
Although the eigenvalues refiect only the noise properties 
and not the signal-to-noise, the growth rate in the fidu- 
cial cosmology is sufficiently flat so that [sA]iA7^'^ shares 
the same form (bottom panel) . Most of the signal comes 
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FIG. 2: Localized constraints on the growth function (j> = Dgiow/a: 
(upper) projected error boxes on the growth function with Icr error 
bars on the Fisher-square root representation with band widths 
taken from the window functions (lower). For comparison the 
growth function for an mtot = 0.2eV model with three equal mass 
neutrinos is shown (dashed line). 



from the first few eigenmodes but the signal-to-noise in 
the first 10-15 eigenmodes remains substantial for a large 
survey. 

Although this principle component analysis is ideal for 
exposing the nature of the information, the oscillating 
windows makes it somewhat difficult to visualize its im- 
pact for model testing. For this purpose, it may be prefer- 
able use more localized linear combinations that retain 
the uncorrelated property at the expense of having over- 
lapping (non-orthogonal) windows. Consider the linear 
combinations defined by rows of the "square root" of the 
Fisher matrix fE2l 



G = SA^/^S* 



(14) 



and with a normalization chosen so that the window el- 
ements sum to unity. As shown in Fig. ^ they yield 
well-localized windows and provide a visualization of the 
data with error boxes whose width is determined as that 
enclosing the central 60% of the window. We show here 
the errors appropriate for a survey with /g^y = 0.1. 

Any model-independent constraints on the growth rate 
can be translated into a model-dependent one by exam- 
ining the of the model fit. In terms of Fisher fore- 
casts, this is equivalent to a re-parameterization through 
the transformation law (^0|). As an example consider a 
family of growth functions that represent a rescaling and 
pivot around a fixed Zp from the fiducial model 



1 + z 
1 + 2:„ 



(f>{id{z) 



(15) 



The pivot point can be chosen to decorrelate the 
errors between A and p by examining Fisher re- 
parameterization of Zp pq|. This choice then has the 



nice property that the error a{A) in the two parame- 
ter model {A,p) are also those of the single parameter 
family of models (A). For the fiducial model and sur- 
vey Zp = 0.64, a{A) = 0.0023/^,^^^ a{p) = 0.0089/r.J^'. 
Physically, such constraints would limit the fraction / of 
the dark matter in a smooth component, for example a 
light neutrino below its free-streaming scale. A smooth 
component induces a change in the growth rate in the 
matter dominated epoch of p = 3/ /5 and a consequent 
change in the amplitude compared with the initial con- 
ditions of 6 A « 4/ (e.g. 0). 

We show a neutrino model with a total mass mtot = 
0.2eV distributed equally into three species Q in Fig. ||. 
This test is potentially substantially more powerful than 
its galaxy clustering analogue due to the lack of an un- 
known bias ^ . While most of the constraint would 
come from the amplitude A, information on the growth 
index p is useful for distinguishing such effects from those 
of the dark energy (§[V) and uncertainties in the initial 
conditions (j|v|). By varying the pixel size one can in 
principle test the scale dependence of growth rate pre- 
dicted in such models. Likewise, deviations would occur 
if the dark energy is not effectively smooth on the scale 
of the pixels. Lensing tomography offers a unique oppor- 
tunity to test the clustering properties of both the dark 
energy and the dark matter on intermediate scales. 



IV. DARK ENERGY DENSITY EVOLUTION 

If the distance-redshift relation D(z) is uncertain due 
to the dark energy, then these degrees of freedom must 
be incorporated into the statistical forecasts as well. For 
simplicity, we will in this context assume that the dark 
matter is composed solely of cold dark matter. 

Fortunately for a smooth dark energy component, both 
the distance-redshift relation (|^) and the growth rate are 
fixed by the dark energy density evolution py)e{z) 
Here the growth rate obeys 
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^ - ^w{z)nDE{z) 
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-[1 - w{z)]n^E{z)^ = , 



(16) 



where the initial conditions are 
and the equation of state 



1 and d(f)/d\iia = 



W{z) ^ ^-O—P, 1- (l*^) 

PDE o a In a 

As in the case of the growth rate, we chose a model- 
independent parameterization as the primary represen- 
tation. Consider the dark energy density in redshift bins 
0, 111, specifically 



di = In 



PcrO 



(18) 
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FIG. 3: Principle components of the dark energy density d = 
lri(PDE/pcro)- (upper) first five eigenfunctions of the growth func- 
tion (thick to thin) plotted in redshift; (lower) rank ordered in- 
verse rms noise X^^^^ (thin) and signal-to- noise ratio in the fiducial 
model. Higher eigenfunctions pick up evolution in the dark energy 
density which is absent in the fiducial model. 



We choose Az of the bins to be the same as the density 
reconstruction Az = 0.1. Since pcro is the critical density 
today, in the limit of fine binning di — In f2DE ■ We also 
require a finite w'{z) = dw/dz to ensure that the dark 
energy remains smooth. Therefore we choose d(z) as a 
spline interpolation of rfj. 

Again, the recovery of finely-binned parameters is 
noisy and correlated across neighboring bins. The infor- 
mation content is best revealed through the Fisher prin- 
ciple component analysis. Shown in Fig. js] (top) are the 
first 5 eigenfunctions. The qualitative difference between 
the eigenfunctions of the growth and that of the density 
is the presence of substantial low redshift information in 
the latter. The dark energy density at low redshift affects 
the distance to all higher redshifts. In Fig. ^ we show the 
eigenvalues as A^^'^ and the signal-to-noise ratio for the 
fiducial model and /sky = 0.1. 

It is again desirable to find an uncorrelated but more 
localized representation of the data. Unfortunately, the 
Fisher square-root technique of Eqn. ( p^ does not yield 
localized windows for the dark energy density. Instead 
we choose a close analogue, the Cholesky decorrelation 
p2t where the windows are the columns of L and Fdd = 
LL* again normalized to sum to unity. The windows are 
shown in Fig. ^ (lower panel). The windows also have 
the interesting property that they are strictly zero below 
some minimum redshift. 

In Fig. m (upper panel) we show the projected con- 
straints on these localized modes compared with the pre- 
dictions for a w = —0.8 model (points) and the actual 
function /Ode(-z) (curve) in this model. The deviation be- 
tween the points and the curve reflects the non-locality 
of the windows. Note that for a cosmological constant 



FIG. 4: Localized constraints on the dark energy density: (upper) 
projected error boxes on the growth function with Icr error bars 
around the fiducial model (line) for the Cholesky representation 
where band widths are taken from the window functions (lower). 
For comparison a model with w = —0.8 is shown both in the predic- 
tions for the modes (points) and the density function itself (curve). 
The small difference reflects the non-locality of the windows. Any 
deviation from constancy represents dark energy that is not a cos- 
mological constant. 



model (pde = const.) there is no deviation by defini- 
tion (straight line) and that the expectation value of all 
points is InfiDE- Any statistically significant difference 
in the values of the points in this reconstruction rep- 
resent a detection of a dark energy component that is 
not a cosmological constant . This is true in spite of 
the non- locality of the windows and independently of the 
model-dependent parameterization of the dark energy. 

Again one may always test specific models for the 
dark energy from the model-independent parameteriza- 
tion. Many dark energy models can be parameterized 
by l^DE, w(zcft ), and w' = dw/dz | z^j,. . The pivot point 
Zcff can be chosen to be the best constrained redshift or 
"sweet spot" by decorrelating the errors in w and w' . As 
was the case for the growth rate, the resulting errors on 
w{zef[) are the same as in the case of a two parameter 
model {iluE,w). For the fiducial model and survey, this 
is Zcs = 0.33 and a{w) = O.OlB/^j^y' . Note that this 
constraint is marginalized over (OdEiW') without prior 
assumptions to their values. 

As in the case of supernovae luminosity distance mea- 
sures ||ll|, there remains a degeneracy between JIde 
and w' in that they may both be adjusted upward to keep 
the dark energy density at the well-constrained low red- 
shifts fixed. In Fig. |5| shows, we show the 68% confidence 
region with various assumptions of prior knowledge on 
r^DE: none, a{nuE) = 0.03, = 0.01. With Ode fixed, the 



errors become cr(w') = 0.046/,i,y' . With no prior, the 
errors degrade to a{w') = 0.069 f^^^' " . 
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FIG. 5: Error ellipses (68% CL) for dark energy parameters (f2DE> 
w, w'): (upper) the (Hde,'^') plane showing the degeneracy direc- 
tion and the efficacy of external information on Hqe; (lower) the 
(OdEi'^) plane with complete and no external information on to'. 
Here w is defined as the equation of state at ^ = 0.33. The fiducial 
model and survey are assumed throughout with /^ity = 0.1. 

Conversely if w' is fixed, the constraints on the dark 

— 1/2 

energy density improve to (T(f7DE) = 0.003/^,^^' . This 
two parameter family (JIde, w) represents the amplitude 
and slope of the dark energy density itself at a normal- 
ization point of z = 0. The remaining degeneracy with w 
can be removed by again going to the "sweet spot" , here 
z = 0.13 where the errors on the dark energy density im- 
prove by a factor of 1.5. Again this also represents the 
errors on fluE in the single parameter family of dark en- 

— 1/2 

ergy models (T(f2DE) — 0.002/gj^y' . Errors on w of course 
remain unchanged. 

V. DISCUSSION 

We have shown that lensing surveys covering more 
than a few percent of the sky with good photometric red- 
shift information can probe the time evolution of the lin- 
ear growth function and distance-redshift relation, both 
of which are sensitive to properties of the dark energy 
and dark matter. Specifically, we have tested a model- 
independent parameterization of the linear growth rate 
and/or dark energy density discretized into bins in red- 
shift. Deviations in the growth rate would indicate a 
component of the dark matter that is not effectively cold 
or dark energy that is not smooth on the lensing scale. 
Deviations in the constancy of the dark energy density 
would rule out a cosmological constant model. 

In this exploratory study, we have made several sim- 
plifying assumptions that would need to be addressed in 
a concrete implementation. Perhaps the primary one is 
that future CMB measurements will completely fix the 
high redshift cosmology. The most uncertain piece in- 



volves the amplitude of the initial fluctuations on the 
scales relevant to the lensing pixels, k ^ 0.05 Mpc~^ 
for degree scales. Fortunately, the pivot point of CMB 
anisotropy experiments with several arcminute scale res- 
olution is sufficiently close to the lensing scale that the 
slope and shape of the initial power spectrum do not 
cause much ambiguity [Q. 

To fully utilize the lensing information, the initial am- 
plitude must be fixed to an accuracy better than the am- 
plitude of the growth function, which we have found to 
be ^ 0.002/^k,y^, i-e. percent level accuracy for surveys of 
several thousand square degrees. For CMB anisotropics, 
this precision requires that the optical depth during 
reionization must be determined to <j{t) 0.01 to resolve 
the amplitude degeneracy. If determined from CMB ob- 
servations alone, this will require polarization measure- 
ments with a precision comparable to the Planck satellite 
( pTf , which can in principle achieve a{hid(^) = 0.0044 
at = 0.05 Mpc~^ 0]. Direct measurements of the 
reionization epoch can also resolve the ambiguity [Q. 
Even in the absence of this information, the evolution of 
the growth and luminosity-distance relation are still con- 
strained. These issues are best addressed through joint 
parameter estimation. 

On the lensing side, the most important assumption is 
that the noise in the convergence measurements is well- 
calibrated and not significantly larger than the projec- 
tions based on intrinsic ellipticities. Even aside from the 
demanding requirements for control of systematic errors, 
there may be intrinsic correlations in the ellipticities 
that need to be modeled or avoided by increasing the 
pixel scale and redshift bin widths. The recovered in- 
formation is largely insensitive to the redshift bin width 
since the high signal-to- noise modes are all low frequency. 

1 /2 

Errors scale roughly as A^'-^ due to the loss of indepen- 
dent modes in a fixed survey area. 

We have also neglected sample covariance between the 
pixels but note that we have correspondingly neglected 
the information contained in such correlations. Indeed, 
we have completely neglected the information contained 
in the non-linear regime which in fact contains the ma- 
jority of the information from lensing tomography [ p^ . 
Clearly, future studies will be required to see how best 
to mine the model-independent information contained in 
lensing tomography. 
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